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space-time dimension 
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A short review of the method for the tensor reduction of Feynman integrals based on recurrence relations 
w.r.t. space-time dimension d is given. A solution of the difference equation w.r.t. d for the n-point one-loop 
integrals with arbitrary momenta and masses is presented. The result is written as multiple hypergeometric series 
depending on ratios of Gram determinants. For the 3-point function a new expression in terms of the Appell 
hypergeometric function Fi is presented. 



1. Introduction 



2. Reduction of tensor integrals 



The increasing precision of current experiments 
and the expected precision of future experiments 
reveal a need to evaluate the two-loop and higher 
order radiative corrections [ An evaluation 
of mass dependent two-loop corrections for many 
physical parameters such as the p parameter [ ^ 
and Z ^bb decay width [ H) is greatly desirable. 
Of special interest is the problem of calculating 
one-loop corrections to the process e+e^ 4 
fermions. To solve this problem efhcient algo- 
rithms for evaluating multi leg one-loop tensor 
integrals are needed. Some progress was recently 
achieved in this direction in [ ^- [ |^. 

In the present paper we consider the algorithms 
based on recurrence relations w.r.t. space-time 
dimension or d-shift recurrence relations which we 
hope will be useful for evaluating multi loop and 
multi leg integrals. In the first part we describe 
the method for the reduction of tensor integrals 
to a set of basic integrals proposed in [ ^ . In the 
second part we present a solution of recurrence 
relations w.r.t. d for one-loop n-point scalar inte- 
grals with arbitrary momenta and masses. This 
is a short report on work done in collaboration 
with J.Fleischer and F.Jegerlehner [||]. 
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The reduction of one-loop tensor integrals to 
scalar integrals with shifted dimensions was given 
in [ 0. The reduction of multi loop tensor inte- 
grals was considered in [ || . In the approach of [ 
p| any tensor L loop integral is represented as 
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ZdL/2 



N 



j=i 

where kj is the momentum of the j-th line, pi are 
external momenta, T^i...a„j^ is a tensor operator 
and 



d 



drrij ' 
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It is assumed that the scalar integral on the r.h.s. 
of (|l|) has arbitrary masses and only after differ- 
entiation w.r.t. raf are these set to their concrete 
values. 

To derive the explicit formula for the operator 
T, L independent auxiliary vectors ai were intro- 
duced to represent the tensor in the integrand as 
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and after that the integral was transformed into 
an a- parametric representation. From the para- 
metric representation it is easy to deduce that 
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where Q is a polynomial in a, a and PiPj. Ap- 
plying T gives a sum made of external momenta, 
g^i, multiplied by scalar integrals with shifted d. 
Advantages of the proposed tensor reduction are: 

• no contractions with external momenta and 
metric tensor and no solution of a linear sys- 
tem of equations are needed 

• it is easy to select a particular tensor struc- 
ture 

• representation in terms of integrals with 
shifted d is compact 

Thus the problem is reduced to the evaluation 
of scalar integrals with different powers of scalar 
propagators and shifts in d. They can be evalu- 
ated by applying generalized recurrence relations 
[ ||. To obtain these relations we start with the 
identity: 
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where R is an arbitrary tensor polynomial. Such 
kind of identities for finding relations between 
Feynman integrals was used in [ pl|,[ [l^ . A sys- 
tematic method was proposed in [ |13[|. Upon dif- 
ferentiating two representations for scalar prod- 
ucts can be used: 

a) integration by parts method [13]. 

kiPi = ^{kf +pf - {ki -pif), 

b) generalized recurrence relations |8| : 
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Combining these methods gives many recurrence 
relations connecting integrals with different in- 
dices Vj and shifts of d. 

2.1. One-loop 7i-point integrals 

The method for the reduction of one-loop ten- 
sor integrals 
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with Cj = {q—pj)'^ — m'j+ie, was presented in [|^. 

The integral In^l can be written as a combination 
of scalar integrals by using the relation 
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Two recurrence relations are needed to reduce 
any scalar integral lli^ to basic ones: 
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where etc. shift the indices i^j Vj ±1 and 
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A„ = A„({pi,mi}, . . .{p„,m„}) 
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Here pi is the external momentum flowing 
through the i-th line and rrij is the mass of the 
propagator corresponding to the j-th line. 

Relations similar to for the integrals with 
zero Gram determinants were given in [ ^ , [ ^ . 

In the next sections we will also use an index 
notation for A„ and G„ 



A„({pii,mjJ, . . . ,{pi„,mi„}), 



A, 



= Gn-liPh,- ■ ■ ,Pi„)- 
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We shall use the index notation for integrals ob- 
tained from I^'^ by contracting some lines. The 
characteristic variables occurring in our deriva- 
tions will be ratios of Gram determinants and 
therefore it is worthwhile to introduce the short- 
hand notation 



fij...k 



Ay.. 



dij ...k 



(8) 



2.2. Two-loop propagator integrals 

The reduction of the two-loop propagator type 
tensor integrals to basic scalar integrals can be 
done by using the relation [ [l^ 



d'^kid'^k^ , ■ , . 
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Ci — k\ — ni\ + ie, C3 — (fci — 



C2 — k 



-it, C4 = {k2 - qf - ml 
C5 = {ki - k2f - ml 



An explicit form for the operator T and all nec- 
essary recurrence relations are given in [ |l^ . The 
recursive procedure allows one to transform any 
diagram into a sum over 30 basic integrals: 
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This algorithm was implemented in [ |l^ as a 
Mathematica package. 

The method of tensor reduction [ ||] was used 
for the evaluation of radiative corrections to sev- 
eral important physical quantities. For example, 
for the evaluation of the 2-loop correction to the 
static potential in QCD [ |l^ and the evaluation 
of the 2-loop correction in the gauge Higgs system 
in 3-dimensions [ ^] . Another important appli- 
cation may be the evaluation of 2-loop correction 
to the Bhabha scattering [ , [ |l§ . 

3. Explicit solution of d shift recurrence 
relations 

The solution of the d— shift relations for one- 
loop propagator type integral was first discussed 
in [|^. Here we present the solution of the relation 
(U) for the integral I^i^ with arbitrary momenta 
and masses and first powers of scalar propagators. 
In this case 
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where the operator k removes the fe-th line from 
If we assume that n — 1 point functions 
are already known then relation ( p^ is an inho- 
mogeneous first order difference equation w.r.t. 
d. Methods of solution of this kind of equations 
are well described in the mathematical literature 
\ M. The redefinition 
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leads to a simpler equation 
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Without loss of generality d can be parametrized 
as 

d=2l- 2e, 

with I being an integer number and e an arbitrary 
small parameter. The solution of equation ( [l^ ) 
then reads 



r=0 



2A. 



Gn-l 



Y^{dMk-li'^-'-'^\ (13) 
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where 6„(e) is an I independent constant. By 
changing bn{e) the solution ( p^ ) can be rewrit- 
ten in another form 



Gn-1 



2A„ 
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The result for 7^''' then reads 
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where 6„ can be determined from the asymptotic 
value of llf' in the limit d — s- cx). 

3.1. Asymptotic form of lli^ at d ^ oo 

The asymptotic form of llf' as d ^ oo can be 
derived from its parametric integral representa- 
tion, lli'^ can be written in a parametric form by 
employing, for example, the following formula 



ClC2 



with 



= r(n) I ... I dxi . . . dxn-i 

■ ■ - Cn Jo Jo 

n-2 n-3 
X^^^2 ■■■^r.-2 ^ 

D = [ciXi. . .X„-l 

-H:2Xi. . .x„_2(l-a::n-i) + - • .+c„(l-a::i)](17) 



then shifting integration variable q and applying 
the formula 



d'^q 
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(18) 



-)" r(a)(mf)"-§ 

Any n-point function will be represented by a 
multiple parametric integral of the form: 

= r f n- " ] ^ • • dxi... dxn-i 



xfi{x}){K{{pj,m,},{x}))^-, (19) 

where {x} = {xi, . . .,Xn-i} and /i„ is a polyno- 
mial. In our opinion this particular parametric 
integral representation is very convenient for per- 
forming asymptotic expansions. The behavior of 
this integral, as d — > oo, can be established by 
using asymptotic methods [ |^ . The main con- 
tribution to the integral will come from the max- 
imum of hn . The maximum can be either on the 
boundary or in the interior of the integration re- 
gion. In case the solution {x} = {x} of the system 
of equations 

dhn 

dxi 

lie in the interior of integration region then it is 
possible to show that 



= 0, 



(1 < i < n - 1) 



(20) 
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and at this point 
hn{,{P],ms},{x}) 



ri...r. 



(21) 



(22) 



In the domain of analyticity ft,„({pj, TOs}, {a;}) > 
and In^'' is an integral of Laplace type. 
In the domain of non-analyticity there are 
subdomains of the integration region where 
hn{{pj,ms}, {x}) < and so the integral J^'^^ has 
an imaginary part. For this kinematic configura- 
tion may be represented as 

{dx}d{h,,)f{{x}) 
{hn{{p„ms}Ax})r-^ 

7T 7T 

[ COS— {d — 2n) + i sin 2" ^ 
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(23) 



5 



Here ftj are subdomains of the integration re- 
gion where /i„ < 0. The boundaries of these sub- 
domains can be found from the solutions of the 
equation 



hn{{pj,ms}, {x}) = 0. 



(24) 



In ( p3| ) we took the main branch of /i„ and as- 
sumed that the momenta and masses are real pa- 
rameters. A concrete example of representation 
( p3| ) will be given in the next section. 
The imaginary part of I^"^ 



n+l) 
{dx}f{{x}) 



\hn{{Pj,ms},{x})\" 



(25) 



is a sum of integrals of Laplace type and therefore 
its asymptotic form at d — > oo can be established. 

It should be noted that representations like (23) 
are valid for multi loop integrals including vac- 
uum ones and also they can be used as an al- 
ternative to the method for evaluating imaginary 
part of the diagram with the help of Cutkosky's 
rules. 

The integrals (|l^) are multiple Laplace type 
integrals of the form 



F{X)^ / f{x)eM>^S{x)]dx, 
Jq 



(26) 



where fi is a bounded simply connected domain in 
k dimensional Euclidean space, x = {xi, . . . Xk), 
A is a parameter, S{x) is a real function. Sup- 
posing that the absolute maximum of S{x) in Q 
is achieved only at the interior point x — x, then 
at A — > cxD 



F(A) - exp[A5(x)] ^ akX^"" 



(27) 



r=0 



This expansion may be differentiated w.r.t. A any 
number of times. The leading term of the expan- 
sion is 



F(A) = exp[A5(x)] (27r/A) 



. f{x) + 0{X-^) 
^y\detS^x{x)\ ' 



(28) 



where |detS'a;a;(x)| is a Hessian determinant. The 
leading contribution from the integral at d ^ oo 



comes from the region where |/i„| has maximal 
value. If an extremum of /i„ occurs inside the in- 
tegration region then we must compare |/i„ | at the 
extremum point ( |2l] ) with its value on the bound- 
ary of the integration region and take the largest 
one. If an extremum of /i„ lies outside the integra- 
tion region then the maximal value of \hn\ will be 
reached on the boundary. Since /i„ in the integral 



is multiplied by the factor 



2 n— 3 



■Xn-2 



the 



non-zero contribution to the leading asymptotic 
comes either from the boundary region = 
or from the boundary region where one of 

Xt = l, {I <i <n-l). 

At these n boundary edges /i„ is equal to hn-i 
of one of the n integrals which one can obtain by 
contracting one of the lines in /„. Similar to /i„ 
the maximum value of these /i„_i is equal either 
to their extremum value 

hn-i = |rii...j„_J (29) 

or to the value on the boundary. Again at the 
boundary will be equal to the hn-2 corre- 

sponding to integrals one can obtain from 
by contracting two lines and so on. Splitting 
the integral for /i„ < into real and imaginary 
parts ( |23|) does not influence our considerations 
because splitting adds a boundary where /i„ van- 
ishes and therefore there will be no contribution 
to the asymptotic value of the integral. 

According to the previous discussion and from 
Eq. (|l^) it follows that in order to find 5„ the con- 
tribution of I^"^ at d ^ oo is needed only if it is 
proportional to r'^^^- Such an asymptotic value 
will occur for kinematic configurations when all 
a; < 1 and /i„ has either maximum or negative 
minimum. In the latter case In^"' must be split 
into real and imaginary parts according to ( [2^ ) 
and the asymptotic value of the term giving a 
contribution proportional to r^^^„ are to be de- 
termined. It should be noted that the contribu- 
tion of order rf^\ may come from the analytic 
continuation of the sums in ( p^ ) . As will be seen 
in the next sections these sums represent hyper- 
geometric functions. 
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From the formula for the leading asymptotic of 
the multiple integral (ESh we obtain: 



-{2n) 



. r 1 



(30) 



and therefore the general solution of Eq. (|1C|) for 
the kinematic configuration satisfying the condi- 
tions < X, < 1 will be 
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dkA, 



^ V 2A,, 

k=l 



E 

a7" 



V 7I'|G„_l 

d - n + 1 
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where (a)^ = T{r + a)/T(a) is the Pochhammer 
symbol. A representation of one-loop integrals in 
terms of Lauricella functions was given in [ . 

4. 2-point function 

Expression (jl^) for llf ^ includes two one-fold 

infinite sums over tadpole integrals (see for- 
mula (ll)): 



(32) 



It is convenient to label the lines of /a'*-' by i,j 
because the integrals 12'*'' will be encountered in 
calculating I^'^'^ , . . . as a result of contraction of 
different lines. 

At n = 2, substituting (|3^ ) into (|5|) gives 



_ '»j \ " \ 2 ), 



^t^-^i^X- (33) 



K) 



One can easily recognize that each of the two in- 
finite sums can be written in terms of Gauss' hy- 
pergeometric function 2F1. The parametric inte- 
gral for /j^' is 



(34) 



with 



I 22 r 2 2 , 2\ , 2 

T'2 = P^jX^ - Xl(p^j - TO, + TO^) + nij. 



An cxtrcmum of h2 is at xi 



TO^)/(2p? ). The maximum of /12 exists if 



' 2p2, < 0, 



(35) 



i.e. only for Euclidean momentum. 

The formula (Eq) gives the imaginary part of 



on the cut 



where 



(d) _ _^ {ix,-x,)ix+~x,)y. 



Pt-' r (^) 



,(36) 



± _ pI -ml^m]± ^-A,j 

is the solution of the equation 

/i2({Py,TOi},a;i) = 0. (37) 

The change of the integration variable x\ = {x\ — 
x'^^z^x'^ makes the integration trivial and gives 



Im/. 



(d) _ 



d-3 



r(d-2) 

Ad) 



Pt 
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The real part of coming from the region 
where /12 < is 

Re/f'=Im4'')ctg|(d-4). 

At large d the asymptotic value of Im /j'''' and 
Re 12'^^ can be easily found. At pfj > (toi + TO2)^ 
both 0<xf <1. 

Now several remarks concerning the evaluation 
of 62 are in order. From (|l5| ) we see that at large 



dimension 62 



d/2 



Such a contribution may 

id) 



come either from the asymptotic value of 
(this may happen only if < < 1) at large d 
or from the analytic continuation of infinite sums 
when their expansion parameter exceeds 1. To 
proceed further let us assume that > m?-. At 
pf - < m'j—mf andp|- > {nii+mjY the value of 62 
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is determined from the asymptotic of /, 



id) _ d/2 

a ■ 

Ad) 



At mj — < < {rrii — rrij) both and 



the infinite sums have no contributions 



d/2 



and therefore in this region 62 = 0. In the region 
{rrii — rtijY < Pij < (m-i + "^2)^ the value of 62 is 
determined from the asymptotic value of infinite 
sums at large d. 

Combining all values of 62 f^nd writing infinite 
sums in terms of hypergeometric functions we ob- 
tained the following result 



r(i-f) 

r(^) - 



di Xij 



d ■ 

2 ' 



d-1 . ^2 



(39) 



The analytic continuation of 2^^! leads to a com- 
pact form for I^'' (see also [ psf) 



i^3((d-2)/2,l,l,l/2,d/2;a;,y) defined as 
F3(a,Q;',/3,/3',7;x,y) = 
^ (a)^(aO»(/3)„.(/30n ^ „ 



7n,n—0 



(42) 



which may be reduced to Fi by means of [ ^ 
Fsia, a' , (3, (3' , a + a' ; x,y) = 

(1 - y)-^'F, (a, f3, f3' , a + a'; x, (43) 
The Appell hypergeometric function ^Fi defined as 

00 

E 



lo i'y)m+nmln\ 



(44) 



Similar to /j'*'* we will label the lines of /g'^^ by 

i, J, A:. The result for /g''^ valid when /13 has max- 
imum inside the integration region reads 

r(d) ' 



r(2 



d\ ^3 



2 = vr y/^g^r^-l 



+Oijk dkXijk + Okij djXijk + Ojki diXijk, (45) 
where 

Xij Oijji — 



5y I- 
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— 2-fl 



r(2-D (^2)1^ 



djXij 



{mfy 



iFi 



1 4-d 

3 . 

2 ' 

r, 2 ' 
3 ■ 

2 ' 
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The expansion of to all orders in e — {A — d)/2 
is given in [ |2^ . 

5. 3-point function 

The expression for l'^^ includes summation 
over 2F1 coming from To derive an explicit 
formula we first use the relation 



(l~z)"2i^l 



— 2^1 



c; z 



1 



,(41) 



in order to remove the factor d from the second 
parameter of 2-^1- Summing over l'^^ gives the 
AppcU hypergeometric functions 



2(rf- 2)^77 

/ 2\—'' 



dpX.j fd-2 dmnj mf 

\ 2 2 2 Tij 1^ Tij 



2{d-2) 
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d-2 ^ I d m^j m'j 
2 2 2 Tij k Tij 



4r(^) 



X2i^l 



1 d~2 . „. . 
2 ' 'ijfc 
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The function Fi has a simple integral represen- 
tation 

/d-2 1 d \ d-2 
^M^'^'2'2'^'^j=^ 



d-4 
U 2 



(1 - xu)^/l - yu 



zdu. 



(47) 
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The polynomial 
/i3 =-xiX2{l -a:i)pi3 -xlx2{l -X2)p\2 

XiX2m\ — Xi[l —Xi){l — X2)P23 

— X2)ni\ + (1 — xi)m\. 
has a maximum at 

91A3 + d2^z diAj 



(48) 



Xi 



-Go 



X2 



diAs + 82 A3 



if 



dx2dx2 



2x\pI^ < 0, 



(49) 



dxidxi dx2dx2 \dx1dx2 

{diA3+d2A3)^ 



2G2 



> 0. 



(50) 



From ([49[) , (|50|) it follows that the maximum will 
be achieved if 



PI2 < 0, and G2 < 0. 



(51) 



The behavior of the integrand for I^'^^ having the 
maximum inside the integration region is demon- 
strated in Fig.l. 




2 . 10 



1 . 10 



Fig. 1: Function xih 



d/2-3 



at 



p?3 = -10, = -5, pI^ = -9, mf = 1, d = 36. 



To evaluate /g'''* in different kinematic regions 
an analytic continuation of 2-F1 and Fi is needed. 
The corresponding formulae are given in [ ^ . The 
representations of I^'^'^ in terms of other hyperge- 



ometric functions can be found in [ E6|-[ 28 



5.1. 4-point function 

The most problematic part of the computation 
of ij^'^^ is the summation of terms from /g'''' with 
Fi function. To simplify the summation we re- 
move the factor d in the first argument of Fi by 
applying the formula [ ^ 

il-xfil-yf'F,ia,P,P\r,x,y)^ 



Fi 7 -a, /?,/?', 7;. 



y 



(52) 



X — 1 y ~ I , 
The result for the 4-point function then reads 



^yki j{d) 



r (2 



where 



3 r(^^) 



+4>ijki diXijki + (jyiijk dkXijki 
+<f>kiij djXijki + 4>jku diXijki, (53) 



T>ijkl 



d-i 
. 2 

ijk 



r,.4; 2F1 



terms with Fi and Fs 



1,^; n,k 

d-2 . 

2 ' 



(54) 



and Fs is the Lauricella-Saran function [ ^ of 
three variables. In our case it has a simple inte- 
gral representation: 



Fs = 



r(f) jy-z)-'^ 



X 



1 arcsin 



l-tz 



dt. 



(55) 



(1 - a; -t- tx)^/l - ty 
All other details one can find in [ |^ . 

6. Conclusions 

The recurrence relations with respect to the 
space-time dimension turn out to be useful for the 
reduction of tensor integrals. The relations for 
master integrals in different dimensions are rather 
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simple even for integrals with several momenta 
and masses. Preliminary investigation of recur- 
rence relations w.r.t. d for the two-loop propa- 14. 
gator type integrals revealed that scalar integral 15. 
in (|^) with all = 1 can be written as multiple 
hypergeometric series. Asymptotic expansion at 16. 
large d may be used as yet another tool for ap- 17. 
proximate evaluation of Feynman integrals. As is 18. 
known such an expansion works well in quantum 
mechanics and in quantum field theory in lattice 19. 
calculations. 

Acknowledgements. I wish to thank orga- 20. 
nizers of the conference "Loops and Legs 2000" 
Johannes Bliimlein and Tord Riemann for the 21. 
useful and well organized conference. I am thank- 
ful to C. Ford and V. Ravindran for carefully 
reading the manuscript and useful remarks. I 22. 
would like also to thank the DFG for financial 
support. 23. 

REFERENCES 

1. J. Erler, S. Heinemeyer, W. HoUik, G. Wei- 

glein and P.M. Zerw as, preprint DESY 00- 25. 
050; |hci"-ph/0005024 

2. R. Barbicri et al. Nucl. Phys. B409 (1993) 

105; 26. 

J. Fleischer, O.V. Tarasovand F. Jegerlehner, 

Phys. Lett. B319 (1993) 249. 27. 

3. J. Fleischer, O.V. Tarasov, F. Jegerlehner and 

P. Raczka, Phys. Lett. B293 (1992) 437. 28. 

4. Z. Bern, L. Dixon and D.A. Kosower, Nucl. 
Phys.B412 (1994) 751. 29. 

5. J.M. CampbeU, E.W.N. Glover and D.J. Mil- 
ler, Nucl. Phys. B 498 (1997) 397. 

6. J. Fleischer, F. Jegerlehner and O.V. Tarasov. 
Nucl. Phys. B566 (2000) 423. 

7. T. Binoth, J.P. Guillct and G. Heinrich, Nucl. 
Phys. B572 (2000) 361. 

8. O.V. Tarasov, Phys. Rev. D54 (1996) 6479. 

9. J. Fleischer, F. Jegerlehner and O.V. Tarasov, 
in preparation. 

10. A.L Davydychev, Phys. Lett. B263 (1991) 
107. 

11. B. Petersson, J. Math. Phys. 6 (1965) 1955. 

12. G. t' Hooft and M. Veltman, Nucl. Phys. 
B192(1981) 159. 

13. F.V. Tkachov, Phys. Lctt.lOOB (1981) 65; 



K.G. Chetyrkin and F.V. Tkachov, Nucl. 

Phys. B192 (1981) 159. 

O.V. Tarasov, Nucl. Phys. B502 (1997) 455. 

R. Mertig and R. Scharf, Comput. Phys. 

Comm. Ill (1998) 265. 

Y. Schroder, Phys. Lett. B447 (1999) 321. 

F. Eberlein, Nucl. Phys. B550 (1999) 303. 

V. Smirnov and O. Veretin, Nucl. Phys. B566 

(2000) 469. 

C. Anastasiou, T. Gehrmann, C. Olcari, E. 



Remiddi and J.B. Tausk, |liep-ph/000326l| . 

L.M. Milne-Thomson, The calculus of Finite 
Differences, Macmillan, London, 1960. 
N. Bleistein and R.A. Handclsman, Asymp- 
totic expansions of integrals Holt,Rinahart 
and Winston, 1975. 

A.L Davydychev, J. Math. Phys. 32(1991) 
1052; 33 (1992) 358. 

F. Berends, A. Davydychev and V. Smirnov, 

Nucl. Phys. B478 (1996) 59. 

A.L Davydychev, Phys. Rev. D61(2000) 087- 

701; A.L Davyd ychev and M.Yu. Kalmykov, 

hep-th/0005287 . 

Higher Transcendental Functions, Vol. 1, A. 
Erdelyi, Ed. (McGraw-Hill Book Company, 
Inc., New York, 1953). 

E.E. Boos and A.L Davydychev, Theor. 
Math. Phys.89 (1991) 1052. 
L.G. Cabral-Rosetti and M.A. Sanchis-Loza- 
no, |hep-ph/9809213| ; 

C. Anastasiou, E.W.N. Glover and C. Oleari, 
Nucl. Phys. B572 (2000) 307. 
S. Saran, Ganita, 5 (1954) 77-91. 



